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Abstract: The accumulation, compression and cooling of the ambient interstellar medium (ISM) in 
large-scale flows powered by OB cluster feedback can drive the production of dense molecular clouds. 
We review the current state of the field, with a strong focus on the explicit modelling and observation 
of the neutral interstellar medium. Magneto-hydrodynamic simulations of colliding ISM flows provide 
a strong theoretical framework in which to view feedback-driven cloud formation, as do models of the 
gravitational fragmentation of expanding shells. Rapid theoretical developments are accompanied by 
growing body of observational work that provides good evidence for the formation of molecular gas 
via stellar feedback - both in the Milky Way and the Large Magellanic Cloud. The importance of 
stellar feedback compared to other major astrophysical drivers of dense gas formation remains to be 
investigated further, and will be an important target for future work. 
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1 Introduction 

The formation of dense, star-forming gas from the dif- 
fuse interstellar medium (ISM) is a key process in the 
evolution of galaxies, and a major unsolved problem 
in astrophysics. While the accumulation, cooling and 
fragmentation of the ISM into self-gravitating, star- 
forming structures has been extensively studied - par- 
ticularly from a theoretical perspective - the subject 
remains a challenging one; involving the interplay be- 
tween a complex array of physical, chemical and astro- 
physical processes. Many key questions remain unan- 
swered, including which astrophysical drivers are re- 
sponsible for the majority of dense gas formation, the 
timescales of cloud formation and destruction, and how 
the details of the fragmentation process set the ba- 
sic foundations for subsequent star formation. In the 
present-day, metal-rich universe, the ISM in star form- 
ing clouds is almost always in the molecular phase. 
The key question can therefore be phrased as: How, 
when and where do we form molecular clouds? 

The basic requirements for molecular cloud forma- 
tion are high densities and high column densities. The 
former is essential for effective cooling (and also fa- 
cilitates rapid molecule formation) and the latter is 
necessary to shield the gas from UV heating (and also 
prevents the photodissociation of molecular species). 
Molecular cloud formation therefore typically involves 
large-scale compressive events capable of piling large 
quantities of material into small volumes. Astrophys- 
ical drivers include global gravitational instabilities in 
galaxy disks (e.g. Wada et al. 2000; Kim et al. 2002; 
Tasker & Tan 2009; Bournaud et al. 2010), the ac- 
cumulation of matter in spiral shocks (e.g. Kim & Os- 
triker 2006; Dobbs et al. 2006; Dobbs & Bonnell 2008), 
and compression in expanding shells driven by stel- 
lar feedback (e.g. McCray & Kafatos 1987; Hartmann 



et al. 2001; Ntormousi et al. 2011); all aided by turbu- 
lence that acts to enhance density on a range of scales 
(e.g. Elmegreen 2002; Mac Low & Klessen 2004; Glover 
& Mac Low 2007; McKee & Ostriker 2007; Federrath 
& Klessen 2012). 

This review will examine the role of large-scale stel- 
lar feedback in the formation and evolution of molec- 
ular clouds. This route to molecular cloud formation 
is of particular interest, since it is key element of how 
star formation self-regulates. The cumulative energy 
input from multiple stellar winds and supernovae may 
form new molecular gas through the sweep up and sus- 
tained compression of the ambient atomic medium in 
giant 'supershells' around OB clusters. Rapid cool- 
ing of this swept-up gas - in combination with grav- 
itational, fluid dynamical and thermal instabilities - 
results in the fragmentation of supershell walls into 
dense clouds, which will become molecular (and even- 
tually star-forming) where density and column density 
requirements are met. This process is often discussed 
in the context of the flow-driven model of molecular 
cloud formation, which envisions dense gas and star 
formation as a rapid, dynamical process, and molecu- 
lar clouds as short-lived entities that form at the in- 
terfaces of turbulent flows (see Walder & Folini 1996; 
Hennebelle & Perault 1999; Elmegreen 2000; Hart- 
mann et al. 2001; Audit & Hennebelle 2005; Vazquez- 
Semadeni et al. 2006, 2007; Hennebelle et al. 2008; 
Heitsch & Hartmann 2008; Banerjee et al. 2009; In- 
oue & Inutsuka 2012; Clark et al. 2012). Magneto- 
hydrodynamical models of this process have been de- 
veloping rapidly in the last decade. 

Despite theoretical advances, the role of large-scale 
stellar feedback in molecular cloud formation remains 
an outstanding question. It is clear that the feedback 
from massive stars plays a key role in both structuring 
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and supporting the disk ISM (e.g. de Avillez & Berry 
2001; Joung et al. 2009; Dobbs et al. 2011; Hill et al. 
2012; Hopkins et al. 2012), and is likely responsible for 
triggering a significant fraction of observed star forma- 
tion on local scales in existing molecular clouds (e.g. 
Boss 1995; Yamaguchi et al. 1999b; Hosokawa & Inut- 
suka 2006; Dale et al. 2007; Leao et al. 2009; Dehar- 
veng et al. 2010). However, its role in the initial pro- 
duction of dense gas is not well constrained. There is 
evidence to suggest that the formation of dense clouds 
and clumps on galactic scales is primarily driven by a 
combination of global gravitational instability and tur- 
bulence, (e.g. Wada & Norman 2001; Elmegreen 2002; 
Tasker & Tan 2009; Bournaud et al. 2010; Mac Low & 
Glover 2012), with stellar feedback relegated to a sec- 
ondary role. In this picture global self gravity produces 
large dense gas complexes and spiral arms, while the 
localised energy input from OB clusters both triggers 
star formation in existing clouds, and drives the for- 
mation of some extra molecular material (Elmegreen 
2002). The magnitude of the stellar feedback contri- 
bution, however, remains to be confirmed. 

Observationally, there is a large body of literature 
that makes at least some reference to the formation, 
evolution and destruction of molecular clouds in re- 
lation to feedback supershells or superbubbles (e.g. 
Jung et al. 1996; Patel et al. 1998; Fukui et al. 1999; 
Kim & Koo 2000; Matsunaga et al. 2001; Yamaguchi 
et al. 2001a,b; Dawson et al. 2008b, 2011a,b). How- 
ever, much of this work consists of case-studies of in- 
dividual objects, in which it is difficult to conclusively 
demonstrate triggered cloud formation, and discussion 
on the origin of the molecular gas is often specula- 
tory. One challenge is that the real ISM is struc- 
turally very complex. It is characterised by a frothy 
structure of loops, filaments, carved out tunnels and 
irregular interlocking shells and bubbles, exhibits frac- 
tal structure across all size scales, and is peppered 
with molecular clouds whose individual evolutionary 
histories are ambiguous. The interaction of super- 
shell shock fronts with pre-existing dense gas is also 
common, and further complicates the interpretation 
of associated clouds. Until recently, the limitations of 
molecular line observations meant that most work on 
the molecular ISM in supershells was carried out in the 
Milky Way, where line-of-sight confusion exacerbates 
observational difficulties. It has therefore proved chal- 
lenging to address the role played by large-scale stellar 
feedback to molecular cloud in a quantitative sense, 
or examine its importance in entire galactic systems. 
Nevertheless, when taken as a whole the literature pro- 
vides a good case for molecular cloud formation in su- 
pershell walls, and there is now a growing body of work 
approaching the problem from a more quantitative or 
statistical standpoint, spurred on by ever-improving 
observational capabilities (e.g. Matsunaga et al. 2001; 
Yamaguchi et al. 2001b; Book et al. 2009; Dawson et al. 
2008b, 2011b). 

This paper reviews the theoretical basis and obser- 
vational evidence for feedback-driven molecular cloud 
formation, focussing primarily on the evolution of the 
neutral ISM and the relationship between supershells 
and the molecular medium. In Section 2 we provide 



some background on the feedback-structured ISM, the 
evolution of a classical shell, and a brief overview of the 
ways in which supershells may influence the molecular 
ISM. Section 3 deals with the theory and modelling 
of molecular cloud formation, including the gravita- 
tional fragmentation of expanding shells, the colliding 
flows model of molecular cloud formation, and numer- 
ical simulations of the ISM in galaxy disks. In section 
4 we touch briefly on the relevant theory pertaining to 
the interaction of a supershell shock front with exist- 
ing dense clouds, although the complex subject of trig- 
gered star formation is given only a brief mention. Sec- 
tion 5 describes observations of molecular gas in feed- 
back superstructures, working outwards from the local 
ISM, to the Milky Way as a whole, to the Magellanic 
Clouds, and presenting strong evidence that molecular 
cloud formation via large-scale stellar feedback is oc- 
curring both in the Galaxy and in the LMC. Finally, 
section 6 provides a brief summary and outlines some 
directions for future development. 

2 Background 

2.1 The Galactic Bubble Bath 

The energy input from massive stars significantly im- 
pacts the structure and evolution of the interstellar 
medium. The ISM of star-forming galaxies is riddled 
with the footprints of this stellar feedback, in the form 
of cool shells, hot bubbles and evacuated cavities seen 
in multiple tracers, from X-ray to radio. The exis- 
tence of such structures in the ISM has been exten- 
sively documented throughout the last century (see re- 
view by Tenorio-Tagle & Bodenheimer 1988), and led 
Brand & Zealey (1975) to coin the phrase 'the cosmic 
bubble bath' to describe the structure of the Galactic 
Disk. The volume filling factor of feedback bubbles in 
the Milky Way is thought to lie somewhere between 
~ 0.05-1.0, with most estimates lying at the lower end 
of this range (Oey & Clarke 1997; Oey et al. 2001; 
McClure-Griffiths et al. 2002; Oey & Garcfa-Segura 
2004; Ehlerova & Palous 2005). 

The term 'supershell' was first used by Heiles (1979) 
to describe atomic hydrogen shells with estimated for- 
mation energies of E > 10 52 erg, but is now generally 
used loosely to refer to any large (R > 100 pc) shell-like 
structure. The theory of the formation of these super- 
structures via correlated supernovae and stellar winds 
was developed during the following years (e.g. McCray 
& Snow 1979; Bruhweiler et al. 1980; Tomisaka et al. 
1981; Mac Low & McCray 1988), and although al- 
ternative formation mechanisms have been proposed 
(e.g. Tenorio-Tagle 1981; Loeb & Perna 1998; Wada 
et al. 2000), stellar feedback is successful in account- 
ing for the observational characteristics, population- 
density and energetics of the majority of objects (e.g. 
McCray & Kafatos 1987; Ehlerova & Palous 1996; McClure- 
Griffiths et al. 2002; Weisz et al. 2009; Warren et al. 
2011). To the present day, a steady stream of obser- 
vational work continues to catalogue and study new 
supershells both in the Milky Way and in external 
galaxies (see e.g. Kim et al. 1999; McClure-Griffiths 



www.pu blish . csiro. au/jo urnals /pasa 



3 




Figure 1: Cartoon showing an edge-on view of the evolution of a supershell in the Galactic Plane (time sequence 
from left to right), illustrating some of the ways in which large-scale stellar feedback can affect the molecular ISM. 
Here black clouds represent molecular gas and the greyscale is the ambient atomic ISM. Labels show examples 
of 1. the triggering star formation in existing molecular gas, 2. the formation of new molecular clouds, 3. the 
disruption and entraining of existing molecular clouds. 



et al. 2002; Ehlerova & Palous 2005; Bagetakos et al. 
2011). 

2.2 Evolution of a Classical Super- 
shell 

The classical analytical model for a stellar wind bub- 
ble expanding into a uniform medium was derived by 
Weaver et al. (1977), and modified for a system formed 
by multiple stellar winds and supernovae by McCray 
& Kafatos (1987). The early evolutionary stages of the 
system are dominated by stellar winds, which create a 
hot, low-density cavity into which subsequent super- 
novae inject their energy. By t ~ 5 x 10 6 yr these 
winds have switched off, and supernovae continue to 
inject energy until t ~ 5 x 10 7 yr, with the mass in the 
cavity acting as a buffer to the supernova blast waves. 
During the stellar wind phase, the idealised bubble 
consists of an inner zone of hypersonic stellar wind, a 
shocked stellar wind layer, a shell of shocked ISM, and 
the ambient ISM. The system initially evolves adiabat- 
ically, expanding much faster than the radiative cool- 
ing timescales in any of the four zones, ending after 
several 10 3 years, when radiative losses in the swept- 
up shell become important and it rapidly begins to 
cool. This drop in temperature is naturally accom- 
panied by a drop in pressure, and the shell collapses 
to a thin, cool layer, compressed by the thermal pres- 
sure of the hot (T ~ 10 6 K) interior. The expansion 
continues to be driven by this thermal pressure, with 
a growth rate of the radius of the shell, R, given by 
R oc i 3//5 . This relation holds until such time as ra- 



diative cooling of the interior gas begins to become 
important, estimated to occur at several 10 6 yr for a 
typical cluster-driven superbubble (McCray & Kafatos 
1987). A shell whose interior thermal energy has been 
entirely radiated away will expand by conservation of 
momentum alone, with R oc t 1 ' 4 . The shell may con- 
tain an ionised inner skin for much of its life, with a 
thickness determined by the radiation output of the 
remaining central sources. 

2.3 Disk Blowout 

Shells with sufficient energy may expand rapidly along 
the Disk vertical density gradient, eventually break- 
ing out entirely and venting their hot interior gas into 
the Halo. Accelerated vertical expansion leads to the 
growth of Rayleigh- Taylor (RT) instabilities, culmi- 
nating in the break-up of the shell and the release of 
its interior gas (Mac Low et al. 1989; Tenorio-Tagle 
et al. 1990). The energy requirements for blowout de- 
pend sensitively on the structure of the local ISM into 
which the central cluster inputs its energy, as a well as 
on the distance of that cluster from the Galactic mid- 
plane, and the strength of the disk magnetic field (e.g. 
Tomisaka 1992, 1998). In a realistic, pre-structured 
medium, however, blowout almost certainly occurs rel- 
atively easily, as an expanding supershell seeks out 
low-density pathways through the inhomogeneous ISM 
(e.g. de Avillez & Berry 2001; de Avillez & Breitschw- 
erdt 2005; Hill et al. 2012). These blown-out 'chim- 
ney' systems form a vital link from the Disk to the 
Halo, supplying the latter with the energy and metal- 
enriched material (Norman & Ikeuchi 1989, see also 
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review by Dickey 2009). 

2.4 The Molecular ISM 

The central question in this review is the role played by 
large-scale stellar feedback on the formation, destruc- 
tion and distribution of the molecular ISM. Figure 1 
shows a cartoon of a supershell expanding in a strati- 
fied medium containing pre-existing molecular clouds, 
which illustrates both the large-scale features of the 
system's evolution, and the range of effects it may 
have on the molecular ISM. An initially atomic shell 
may persist for long enough and accumulate sufficient 
material to become dense, cold and molecular, eventu- 
ally fragmenting to form discrete molecular clouds and 
new stars. A shell may also encounter existing dense 
clouds in its path, resulting in their dynamical disrup- 
tion and eventual destruction. Such encounters are 
also quite capable of triggering star formation in pre- 
existing dense clouds, however, so that the formation 
and subsequent collapse of fresh molecular material is 
not the only route to star formation afforded by large- 
scale stellar feedback. While we do not focus on the 
question of triggered star formation in this review, it is 
prudent to point out that the 'destruction' of a molec- 
ular cloud can not automatically be assumed to have 
a negative impact on star formation. 

3 Molecular Cloud Formation 
in Supershells: Theory & Mod 
elling 

An extensive body of literature exists on the physics 
and chemistry of the neutral ISM, and on the transi- 
tion between its atomic and molecular phases. Here we 
focus on that work which is of particular relevance to 
the formation of molecular clouds by large-scale stellar 
feedback - either through direct modelling of super- 
shell systems, or through the development of general 
theory of ISM flows that includes those driven by su- 
pernovae and stellar winds. 

3.1 Molecule Formation & Destruc- 
tion 

While recent work argues that the presence of molecules 
is not necessary for the formation of cold, dense gas 
(Glover & Clark 2012; Mac Low & Glover 2012), molec- 
ular emission lines are the primary tracer of the star- 
forming ISM, and an understanding of how the major 
species form and evolve is essential in interpreting ob- 
servational data. From a theoretical perspective, a re- 
gion of the ISM in which hydrogen is predominantly in 
the form of H2 constitutes a molecular cloud. However, 
observationally it is the abundance and properties of 
trace molecules - in particular CO - that determine 
whether a cloud is detectable. 

Most treatments that explicitly follow molecule for- 
mation make use of approximate expressions for the 
formation rate of H2 on dust grains (e.g. Hollenbach & 



Salpeter 1971; Tielens & Hollenbach 1985). For condi- 
tions typical of the cold neutral medium (CNM, T ~ 
80 K, n ~ 50 cm" 3 , see e.g. Field et al. 1969; Ferriere 
2001), the characteristic formation rate is ~ 3 x 10 7 
yr, but a strong density dependence ensures that H2 
forms rapidly for densities that are even moderately 
enhanced with respect to canonical CNM values (see 
also Glover & Mac Low 2011). For H2 to survive, it 
must be shielded from dissociating photons in the en- 
ergy range 11.08 < hv < 13.6 eV. While full modelling 
of H2 dissociation and shielding is complex (see e.g. 
van Dishoeck & Black 1988; Draine & Bertoldi 1996), 
a combination of efficient self-shielding and dust shield- 
ing means that for UV field strengths within a factor 
of few of the Draine field, hydrogen will typically be in 
its molecular form for visual extinctions of Av > 0.1- 
0.5 (e.g. Franco & Cox 1986; van Dishoeck & Black 
1988; Bergin et al. 2004; Wolfire et al. 2010). This 
corresponds to column densities of roughly ~lx 10 21 
cm" 2 (Draine & Bertoldi 1996). 

For CO, the most commonly used line tracer of the 
molecular ISM, formation proceeds efficiently through 
collisional processes (van Dishoeck & Black 1988), but 
less efficient self-shielding means that abundances are 
primarily determined by visual extinction (Glover & 
Mac Low 2011). The shielding requirements are there- 
fore somewhat more stringent than for H2, at Av > 
0.6-1.0 (van Dishoeck & Black 1988; Bergin et al. 2004; 
Wolfire et al. 2010). An interesting implication of 
this is the existence of a substantial quantity of 'dark' 
molecular gas not seen in the usual surveys (e.g. Gre- 
.nier et al. 2005; Wolfire et al. 2010). 

3.2 Gravitational Instability of Ex- 
panding Shells 

For molecular clouds to fulfil their role as star forma- 
tion sites, they must become self-gravitating - or at 
least contain self-gravitating substructure. Motivated 
by this, a number of authors have applied the theory 
of gravitational instabilities in expanding shells to the 
parameter space appropriate to supershells. McCray & 
Kafatos (1987) use an analytical approximation for the 
fragmentation of an expanding shell into gravitation- 
ally bound clouds to derive a time for the onset of grav- 
itational instabilities as t « 3.2 x 10 7 N* n 1 ^ 2 a 5 J S 
yr, where N* is the number of supernova progenitors 
in the parent cluster, no is the density of the ambi- 
ent medium and as is the magneto-sonic speed in the 
shell in km s _1 . For typical supershell parameters this 
suggests a timescale of a few 10 7 years for shells to 
fragment into gravitationally bound clouds. The au- 
thors also state that they expect supershells to be- 
come predominantly molecular well before this, within 
timescales of ~ 10 6 yr. However, this appears to be 
far too short for shells to accumulate sufficient mate- 
rial for H2 shielding, and is in general not consistent 
with other models or with observations. 

A series of numerical models by Ehlerova et al. 
(1997), Efremov et al. (1999), Ehlerova & Palous (2002) 
and Elmegreen et al. (2002) use the thin shell approx- 
imation (Sedov 1959; Kompaneets 1960) to study the 
gravitational fragmentation of supershells over a broad 
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range of parameter space, with particular focus on 
galactic environment. It should be noted that recent 
work by Dale et al. (2009) and Wiinsch et al. (2010) 
has demonstrated that the thin shell approximation 
deviates from the predictions of numerical simulations 
in cases where the external confining pressure is small. 
Nevertheless, it provides a useful and computationally 
cheap method of exploring the large regions of param- 
eter space appropriate for supershells evolving under a 
wide range of different conditions. 

Ehlerova et al. (1997) note that fragmentation timescales 
are insensitive to the initial energy input, but are strongly 
density dependent, and assume a realistic input clus- 
ter size of 40-100 supernova progenitors to derive the 
condition that systems with n < 0.3 cm~ 3 never be- 
come unstable. Similarly, Ehlerova & Palous (2002) 
derive a critical column density for gravitational frag- 
mentation of TV > 10 20 -10 21 cm for realistic values 
of the energy input and sound speed in the ambient 
medium. Characteristic fragmentation timescales for 
solar neighbourhood densities of no ~ 1 cm~ 3 are long, 
at ~ 5 x 10 7 yr (Ehlerova et al. 1997), but would be 
significantly faster for moderately over-dense regions. 
Similarly, in a stratified medium, it is the dense central 
regions within the galactic plane that become unsta- 
ble to gravitational collapse, while polar regions often 
remain stable throughout a shell's lifetime (see also 
Mashchenko & Silich 1994). Here, the thickness of the 
disk is critical in determining the fate of the system, 
with Gaussian scaleheights of a < 100 pc resulting 
in shells that never fragment. Galactic rotation also 
strongly affects the evolution of the shell. Shear de- 
forms shells into elongated ellipses, and mass accumu- 
lated in the shell slides to the tips, forming instability 
regions there. This phenomenon has also been noted 
by Tenorio-Tagle & Palous (1987), who demonstrate 
that in a typical Milky Way environment these sites 
can accumulate sufficient matter to go molecular and 
form giant molecular clouds (GMCs). The scenarios 
explored in these studies imply that the type of galaxy 
into which a supershell is born strongly influence its 
propensity for gravitational instability, and hence the 
likelihood of forming new molecular clouds through 
this mechanism. Elmegreen et al. (2002) derive a set of 
dimensionless conditions for this 'triggering', and find 
that dwarf galaxies such as the LMC have an advan- 
tage over early type spirals such as the Milky Way. 

3.3 Molecular Cloud Formation in 
Colliding Flows 

A realistic, turbulent ISM is subject to a range of pro- 
cesses that are not treated in the simplified gravita- 
tional stability analysis above, yet which may have a 
profound effect on its evolution. The last 10-15 years 
have seen a profusion of hydrodynamic and magneto- 
hydrodynamic (MHD) simulations that attempt to model 
the ISM in increasingly realistic ways, many focussing 
on the formation of molecular clouds. 

Much work centres around the paradigm of flow- 
driven molecular cloud formation, in which the produc- 
tion of star- forming molecular gas is integrated into the 
modern picture of a dynamic, turbulent ISM. In this 



picture, the compression, cooling and fragmentation 
of the atomic medium in colliding ISM flows produces 
clumpy sheets and filaments of cold, turbulent mate- 
rial, which go on to become highly-structured, self- 
gravitating molecular clouds once density and column 
density requirements are met (see review by Vazquez- 
Semadeni 2010). Hartmann et al. (2001) lay out part 
of the astrophysical motivation for this picture, which 
has its roots in observational evidence that molecu- 
lar clouds must form rapidly, birth stars rapidly, and 
disperse rapidly (see also Elmegreen 2000; Hartmann 
2003). These authors use simple analytical approxi- 
mations to argue that material accumulated in large- 
scale ISM flows will become molecular, magnetically 
supercritical and self-gravitating on roughly similar 
timescales, and that supershells formed by stellar feed- 
back are excellent candidates for the drivers of such 
flows. 

A combination of thermal and hydrodynamical in- 
stabilities, together with turbulent compression (see 
e.g. Federrath et al. 2010), is critical in condensing 
cold gas from the warm ambient medium. The atomic 
ISM is a thermally bistable medium (Field et al. 1969; 
Wolfire et al. 1995), with two linearly stable phases 
corresponding to the warm neutral medium (WNM) 
and cold neutral medium (CNM) - the latter being the 
atomic precursor to molecular clouds. Trans-sonic con- 
verging flows of initially stable WNM trigger runaway 
cooling and the formation of cold gas (Hennebelle & 
Perault 1999), with dense structures developing rapidly 
on sub-parsec scales (e.g. Heitsch et al. 2006; Vazquez- 
Semadeni et al. 2006). Models have progressed to in- 
clude increasingly realistic physics and explore various 
different aspects of the cloud formation process, in- 
cluding the roles of turbulence (Audit & Hennebelle 
2005; Vazquez-Semadeni et al. 2006; Glover & Mac 
Low 2007), magnetic fields (Hennebelle et al. 2008; 
Inoue & Inutsuka 2008; Banerjee et al. 2009; Heitsch 
et al. 2009; Inoue & Inutsuka 2009), molecular chem- 
istry (Bergin et al. 2004; Glover & Mac Low 2011; 
Clark et al. 2012; Glover & Clark 2012; Inoue & Inut- 
suka 2012) self-gravity (Vazquez-Semadeni et al. 2007; 
Heitsch et al. 2008b,a; Hennebelle et al. 2008), and the 
interplay between different instabilities (Heitsch et al. 
2008a). 

Taken together, these simulations suggest that for 
expansion velocities and ambient densities typical of 
Galactic supershells (v exp ~ 10-20 km s" 1 and no ~ 1- 
5 cm -3 ), substantial quantities of CO-rich molecular 
gas can be produced on timescales of a few 10 6 to 
~ 10 7 yr (Bergin et al. 2004; Heitsch & Hartmann 
2008; Clark et al. 2012; Inoue & Inutsuka 2012), pro- 
vided that magnetic fields do not provide significant 
support against collapse (Inoue & Inutsuka 2008, 2009; 
Heitsch et al. 2009). For the parameter space appro- 
priate to the atomic ISM, cooling timescales are much 
shorter than the characteristic timescales for gravi- 
tational instability, and thermal/dynamical instabili- 
ties dominate over gravity in driving fragmentation in 
the compressed medium (Heitsch et al. 2008a). This 
implies that small, dense, thermally-driven condensa- 
tions will develop long before a supershell becomes 
gravitationally unstable, and that some molecular gas 
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Figure 2: High resolution, two-dimensional hydrodynamical simulations of cold gas formation at the interface 
of two superbubbles colliding in a turbulent diffuse medium (reproduced with permission from Ntormousi et al. 
2011). These snapshots show the evolution of the system at 3 Myr (left) and 7 Myr (right) after the start of the 
simulations. The top panels show hydrogen number density the bottom panels show gas temperature. After 7 
Myr, copious amounts of clumpy and filamentary cold gas has been formed. 



can form without the ISM becoming self-gravitating 
(Heitsch & Hartmann 2008). However, global contrac- 
tion helps to reduce overall timescales, and in partic- 
ular may be important in accumulating sufficient ma- 
terial to shield CO so that clouds become observable 
(Heitsch & Hartmann 2008). Once initiated, gravita- 
tional collapse will proceed hierarchically, beginning 
first in the dense substructures already imprinted on 
the medium by thermal instabilities and turbulence 
(Heitsch et al. 2006; Vazquez-Semadeni et al. 2007). 
One implication of this is that clouds are expected to 
begin forming stars rapidly once they become visible 
in molecular tracers, suggesting that the triggered for- 
mation of molecular clouds is synonymous with the 
triggered formation of stars. 

These results provide a strong theoretical basis for 
molecular cloud formation in supershells, and demon- 
strate its viability in systems with flow properties that 
are consistent with observed objects. The next stage is 



to explicitly simulate supershell systems as opposed to 
generalised flows. So far only one high-resolution hy- 
drodynamical simulation has explicitly modelled flow- 
driven molecular cloud formation in a supershell sys- 
tem. Ntormousi et al. (2011) present 2D models of two 
superbubbles expanding into a uniform homogeneous 
or turbulent medium, and investigate the formation of 
dense gas around the bubble peripheries and interac- 
tion zone. Their supershells are blown by supernovae 
and time-dependent stellar winds with properties cal- 
culated from population synthesis models, and evolve 
with realistic flow velocities, timescales and size scales. 
Figure 2 reproduces their snapshot of the temperature 
and density of the turbulent model at times of 3 and 7 
Myr. By the end of a 7 Myr run, their simulation box is 
filled almost entirely by the shell systems, and a com- 
bination of non-linear thin shell, thermal and Kelvin- 
Helmholz (KH) instabilities has lead to the formation 
of copious amounts of clumpy and filiamentary cold 
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(T < 100 K) gas - mostly associated with the shell 
collision zone. They find that ~ 65-85% of the gas is 
contained in these small, dynamical structures, which 
have characteristic sizes of < 1 pc and densities of 
10 2 -10 3 cm - fulfilling the approximate density and 
column density requirements for molecule formation. 

The two-dimensional nature of these simulations 
imposes some limitations on the evolution and struc- 
ture of the turbulence and fragmentation. Unlike in 
3D models, where large-scale structures break up into 
smaller fragments, in 2D simulations they have a ten- 
dency to merge into larger agglomerations (e.g. Fed- 
errath et al. 2010), potentially affecting the ease with 
which dense clouds are formed. Conversely, the lack 
of self-gravity in the simulations is expected to have 
the opposite effect - making it harder to form dense 
clouds. The inclusion of magnetic fields in future mod- 
els would also modify the properties of the collision 
zone (e.g. de Avillez & Breitschwerdt 2005) and likely 
provide additional support. Nevertheless, these mod- 
els are a promising first step, and will be developed 
extensively over the coming years to include more re- 
alistic physics. 

Finally, it is likely that both pre-existing inhomo- 
geneities in the ISM and magnetic field orientation lead 
to a selection effect for molecular cloud formation in 
supershell walls. A moderate elevation in the mean 
density along a swept-up column significantly reduces 
cooling and molecule formation timescales, suggesting 
that cloud formation sites may be determined by the 
placement of pre-existing concentrations of moderately 
over-dense gas such as CNM sheets or filaments (see 
also Dobbs et al. 2012, for similar findings on galactic 
scales). Conversely, magnetic pressure has the power 
to completely prevent the formation of molecular gas 
when a significant component of the field exists per- 
pendicular to the flow direction, leading to the re- 
quirement that this perpendicular component be van- 
ishingly small (Inoue & Inutsuka 2009; Heitsch et al. 
2009). While it remains to be seen how the inclusion of 
self-gravity and non-ideal magneto-hydrodynamic pro- 
cesses affect the rigidity of this conclusion, the selec- 
tion of cloud formation sites by field orientation and 
mean density may provide an attractive explanation 
for why real supershells are not more molecular than 
they are. 



ken fragments with 'tunnels' channeling their hot inte- 
rior gas out into the halo (e.g. de Avillez & Berry 2001; 
de Avillez & Breitschwerdt 2005; Joung & Mac Low 
2006; Joung et al. 2009; Hill et al. 2012). An example 
is shown in Figure 3, from Hill et al. (2012). While 
these models lack the resolution to properly resolve 
very dense gas formation, the fact that the bulk of the 
cold material is concentrated into high density clouds 
in shock-compressed layers around or between bubbles 
(e.g. de Avillez & Breitschwerdt 2005) is highly con- 
sistent with the scenario described in §3.3. 

The lack of self-gravity and realistic disk dynam- 
ics mean that while the above models are a promising 
demonstration of the feasibility of molecular cloud for- 
mation in stellar feedback flows, they are unable to 
address the critical question of its relative importance 
in a galactic context. Simulations of rotating, self- 
gravitating galaxy disks provide a more complex pic- 
ture, in which large-scale gravitational instabilities and 
spiral arms play a major role in dense gas formation, 
and the contribution of stellar feedback is by no means 
clear (e.g. Wada & Norman 2001; Shetty & Ostriker 
2008; Bournaud et al. 2010; Tasker 2011; Dobbs et al. 
2011, 2012). Recently, Dobbs et al. (2012) have explic- 
itly addressed the question of which mechanisms are 
responsible for driving the compressive motions that 
form the GMCs in their simulations, and find that 
the primary drivers depend on the properties of their 
model galaxies, such as the nature of the spiral poten- 
tial and the level of star-formation feedback. While 
conclusions remain open, a future in which these pre- 
dictions can be directly compared with observations is 
promising, as new facilities such as ALMA come on- 
line. 



4 Pre-Existing Molecular Clouds 

Supernovae and stellar winds propagate into a highly 
structured, inhomogeneous medium that already con- 
tains dense gas. The effects of cluster feedback on ex- 
isting molecular clouds range from the violent ionisa- 
tion, compression and dissipation of the parent cloud 
at early evolutionary stages, to the weak interaction 
of a cloud with an old, almost-stalled supershell. The 



3.4 



question of how feedback affects parent molecular clouds 
Whole-Disk Models of the Feedback-is a complex one, involving the combined effects of 



Structured ISM 



While high-resolution numerical simulations have ex- 
amined in detail the small-scale processes associated 
with gas cooling and fragmentation, another class of 
model tackles the global effects of stellar feedback on 
the structure and dynamics of the ISM. Such simula- 
tions generally model kpc-scale sections of galaxy disk 
at moderate resolution with realistic supernova input 
rates, and have been very successful in producing the 
cold, warm neutral, warm ionized, and hot components 
of the ISM with reasonable mass fractions, distribu- 
tions and structure. Salient features include a thin cold 
disk, a 'frothy' disk of warmer material, and the exis- 
tence of a great many shells, supershells and their bro- 



strong ionisation, radiation pressure and mechanical 
energy input on a pre-structured cloud. The interplay 
between these processes and the importance (or not) of 
triggered star formation is still a topic of much debate 
and will not be dealt with in this review (see e.g. Dale 
et al. 2012; Walch et al. 2012, and references within). 
Instead we focus primarily on the region of parame- 
ter space appropriate to a system that has expanded 
far beyond its birth site and out into the surrounding 
medium - i.e. an object that would already be defined 
as a 'supershell' by an observer. For an idealised shell 
of R > 50 pc this implies expansion velocities of no 
greater than a few tens of km s _1 (McCray & Kafatos 
1987); translating to shocks that are mildly supersonic 
in the WNM. 
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Figure 3: Images of density and temperature from magnetohydrodynamic simulations in which supernovae drive 
turbulence and establish a multi-phase, stratified medium (reproduced with permission from the erratum to Hill 
et al, 2012). The images are two-dimensional slices through a three-dimensional simulation, where the x axis 
is in the midplane and the z axis shows the vertical distance from the midplane. The model used here is the 
magnetised model with 2 pc resolution in the midplane ("bx50hr") described by Hill et al. (2012). Energy injected 
by supernovae creates hot, low density remnants; surrounded by dense, long-lived filaments of cold gas. 



4.1 Cloud Disruption 

The dynamical interaction of a dense cloud with a 
shocked flow is commonly parameterised in terms of 
the 'cloud crushing time', t cc = (ro/vi)(n c i/m)° , where 
ro is the cloud radius, Vi is the velocity of the shock 
in the ambient medium, and n c i and i%i are the num- 
ber densities of the cloud and the ambient medium re- 
spectively. The simplest case is a steady planar strong 
shock and a smooth spherical cloud, in which radiative 
losses, thermal conduction, gravity and magnetic fields 
are neglected. Such clouds are significantly disrupted 
within a few t cc , on a destruction timescale tdest, de- 
fined as the time taken for the core mass of the cloud to 
reduce be a factor of 1/e (Klein et al. 1994; Nakamura 
et al. 2006). For the test case of a small (ro = 2 pc), 
moderately dense (n c ; = 200 cm -3 ) cloud encounter- 
ing a gently expanding (vi = 20 km s _1 ), thick, dense 
(m = 10) shell, t cc « 0.4 Myr, implying t de st < 2-3 
Myr. 

In reality this is likely a lower limit to cloud sur- 
vival times. The Mach numbers of supershell shocks 
are relatively low {M ~ 2.5 for a 20 km s _1 shock 
in the WNM), and timescales may be several times 
longer than the strong shock case (Nakamura et al. 
2006; Pittard et al. 2010). Radiative cooling also can- 
not be ignored for supershell-molecular cloud interac- 
tions. Its inclusion inhibits cloud destruction, encour- 
ages the formation of over-dense clumps and filaments 
and can significantly extend timescales above the adi- 
abatic limit (Mellema et al. 2002; Fragile et al. 2004; 
Orlando et al. 2005). The role of magnetic fields is 



more complex - while they tend to inhibit the develop- 
ment of hydrodynamic instabilities and reduce mixing, 
both orientation and field strength are important, and 
some configurations may result in more efficient frag- 
mentation of the cloud material (Gregori et al. 1999; 
Orlando et al. 2008; Shin et al. 2008). Conversely, en- 
vironmental turbulence has been shown to speed up 
cloud destruction, though the magnitude of the effect 
is sensitive to the properties of the assumed turbulence 
and is still under investigation (Pittard et al. 2009, 
2010). Nevertheless, it seems reasonable to assume 
that dynamical disruption of molecular clouds by an 
evolved shell proceeds reasonably slowly, on timescales 
comparable to the shell lifetime. It is worth noting, 
however, that the physical stripping and dissipation of 
cloud material also renders a disrupted cloud increas- 
ingly vulnerable to the UV dissociation of CO. The 
observable lifetime of the entity we see as a CO cloud 
may therefore be shorter than the survival time of the 
dense material itself (see also Dawson et al. 2011b). 

The efficiency of momentum transfer is also impor- 
tant, because it determines how readily the shocked 
cloud can be accelerated by the expanding shell. Adi- 
abatic models find that characteristic drag timescales 
are generally comparable or lower than tdest, suggest- 
ing that a shocked cloud - or its remaining fragments - 
will undergo significant acceleration before being com- 
pletely disrupted (Klein et al. 1994; Nakamura et al. 
2006). Higher density contrasts and lower Mach num- 
bers result in less efficient acceleration (e.g. Pittard 
et al. 2010), as does any process that reduces the sur- 
face area of the cloud perpendicular to the shock di- 
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rection (e.g. radiative cooling; Orlando et al. 2008), 
while environmental turbulence and certain magnetic 
field orientations may increase it (Shin et al. 2008; Pit- 
tard et al. 2009). As in the case of cloud destruction, 
the details of the cloud acceleration depend both on 
the model physics and on the cloud-shock parameters, 
and the final velocity of the accelerated cloud material 
may be anywhere between ~ 0.1-0.8 times the initial 
shock velocity. 

The majority of these studies assume the post- 
shock flow has an infinite depth - a situation that is 
clearly inappropriate for the shell-cloud case. The fi- 
nite depth of a swept-up shell limits the time of the 
flow-cloud interaction, and a cloud encountering a su- 
pershell will likely pass into the hot shell interior be- 
fore it is destroyed, leaving a tail of shell material 
trailing behind it (Pittard 2011). For a shell thick- 
ness of ~ 10 pc, the example cloud described above 
will enter the interior on timescales of ~ 1 to several 
Myr, depending on the drag efficiency. In the interior 
regime the material flowing past it will be more diffuse 
and less dynamically disruptive, but will likely be hot 
ionised medium. The classical evaporation timescale 
for clouds embedded in fully ionised medium is given 
by t evap ~ 3.3 x 10 20 n cl Tr 5/2 r 2 cliPC yr (Cowie & Mc- 
Kee 1977), which equates to ~ 10 8 yr for the example 
cloud, assuming Ti = 10 6 K. This suggests that ther- 
mal evaporation is unlikely to be a dominant destruc- 
tive influence. 

Finally, while we do not deal explicitly in this re- 
view with the topic of triggered star formation, it should 
be stressed that the destruction of a molecular cloud 
can not be assumed to have a negative impact on star 
formation. Indeed, for moderate shock velocities ap- 
propriate to supershell systems, star formation may be 
triggered readily for a range of molecular cloud proper- 
ties (e.g. Vanhala & Cameron 1998; Melioli et al. 2006; 
Leao et al. 2009). These considerations are important 
in interpreting observations of supershell-associated molec- 
ular gas and stars. 

5 Observations of Molecular Gas 
in Supershells 

5.1 General Considerations 

In this section we will review observational evidence 
for molecular cloud formation in supershells, as well 
as some broader examples of the interaction between 
stellar feedback and the molecular ISM. 

The most commonly-used spectral line tracers of 
the molecular ISM are the low- J transitions of CO. CO 
is the workhorse of molecular ISM astronomy, since its 
high abundance (nco/nn 2 ~ 10~ 4 ), low-lying rota- 
tional energy levels (AE/k ~ 5 K for J=l) and rel- 
atively low critical densities (n ~ 1000 cm~ 3 for the 
lower J transitions) mean that it is readily observable 
even in relatively diffuse, quiescent molecular gas. For 
a detailed physical and chemical census of a molecular 
cloud it is possible to assemble a suite of lines from 
multiple tracers that probe a wide range of density 



and temperature regimes; an approach that is partic- 
ularly useful for detailed studies of the star formation 
process. However, from a molecular cloud formation 
perspective we are generally more concerned with iden- 
tifying zones in which the ISM is simply in its molec- 
ular form. CO is an excellent choice for this, although 
it is worth noting that even CO fails to trace the most 
diffuse molecular gas, where hydrogen is in the form 
of H2 but carbon is atomic (e.g. Reach et al. 1994; 
Grenier et al. 2005; Wolfire et al. 2010). 

Molecular gas occupies only a small volume frac- 
tion of the ISM, and probes of other physical regimes 
are needed to form a complete picture of a supershell 
system (see e.g. Heiles et al. 1999) . Hi observations in 
particular are indispensable in probing the structure 
and kinematics of the neutral ISM, and in correctly 
relating molecular clouds to the global structure of a 
supershell system. Much of the work discussed below 
combines Hi and CO data to investigate the structure 
and evolution of the ISM in supershells. Other tracers 
such as Ha and soft X-rays are also useful, and provide 
valuable information on the ionised inner walls and the 
hot diffuse medium in shell interiors. 

A key morphological/kinematic indicator of poten- 
tial cloud formation is co-moving CO clouds that form 
coherent parts of an atomic shell - either lying along Hi 
walls or well embedded within them. The converse is 
molecular clouds that show a head-tail morphology or 
other signs of dynamical disruption, or are entrained in 
a shell interior, likely indicating pre-existing dense gas. 
This simple diagnostic is by no means perfect. Various 
hybrid scenarios in which molecular clouds are formed 
early in a system's evolution and later impacted by 
subsequent supernova blast waves are certainly possi- 
ble. Moreover, analysis can be frustrated by the ir- 
regularity of real observed supershells and line of sight 
confusion. Nevertheless, morphological diagnostics are 
an excellent starting point. 

A complementary approach is to measure whether 
supershells are associated with a net increase in the 
molecular fraction of the ISM in the volumes they oc- 
cupy. This attempts to answer the question of whether 
the formation of new molecular gas has dominated over 
the destruction of pre-existing clouds, as measured at 
a particular epoch in the system's evolution. This is 
an approach taken by Dawson et al. (2011b, 2012b, 
discussed below) by comparing the H2 mass fraction 
in supershell volumes with nearby control regions. 

Finally, stellar population studies are a powerful 
tool for reconstructing the history of a system. As 
well as providing valuable constraints on the original 
energy input from the central cluster, they can also 
provide insight into the nature and timing of any trig- 
gering that may have occurred during a supershell's 
evolution. 

5.2 The Local ISM 

The local ISM, within a few hundred parsecs of the 
Solar System, is in many ways a microcosm of much 
of the Milky Way Disk. It is an environment that has 
been carved out and shaped by multiple generations 
of stellar activity - characterised by a complex struc- 



10 



Publications of the Astronomical Society of Australia 



ture of loops, filaments, tunnels, irregular shells and 
bubbles; and populated by the stellar clusters and as- 
sociations variously responsible for and triggered by 
this activity (see de Geus 1992; Heiles 1998; Lallement 
et al. 2003; Frisch et al. 2011). Unsurprisingly, the dis- 
tribution of the local molecular ISM is closely related 
to these superstructures. 

Perhaps the most striking example is the Gould 
Belt and its Hi counterpart the Lindblad Ring. The 
Gould Belt /Lindblad Ring is the expanding, inclined 
ring of neutral gas and star-forming regions that con- 
tains all of the major local molecular cloud complexes 
within a distance of ~ 500 pc from the Sun, including 
Orion, Taurus, Perseus, Ophiuchus and Lupus (Tay- 
lor et al. 1987; see also review by Poppel 1997). If the 
Gould Belt is a feedback structure, it is a superb candi- 
date for the formation of molecular gas and stars from 
a swept-up supershell. Indeed, one scenario holds that 
the Gould Belt clouds are the remnants of an old fossil 
supershell that has fragmented sometime in the last 
15-25 Myr to form gravitationally bound star-forming 
complexes (Olano 1982; see also Bally 2001). However, 
while stellar feedback is a strong contender, the origin 
of the Gould Belt is unclear, with other candidates in- 
cluding the impact of a high-velocity cloud (Comeron 
& Torra 1992, 1994) and the braking of a supercloud 
entering the spiral arm (Olano 2001; see also Grenier 
2004 and references within). 

As this uncertainty suggests, reconstructing a re- 
gion's history is often a challenging exercise in astro- 
nomical forensics. Nevertheless, there are numerous 
indications that stellar feedback has played an impor- 
tant role in the evolution of the local molecular ISM. 
We will explore a few illustrative examples in the para- 
graphs to follow. 

The Lupus and Ophiuchus molecular cloud com- 
plexes form part of the neutral ISM delimiting the cav- 
ity of the Local Bubble - the old supershell in which 
the Sun is currently located (see Lallement et al. 2003). 
Their distances are ~ 155 pc and ~ 120 pc respectively 
(Lombardi et al. 2008), and both complexes have esti- 
mated (CO-bright) molecular masses of ~ 10 4 M (de 
Geus & Burton 1991; Tachihara et al. 2001), meaning 
that they are smaller than typical GMCs. The Lu- 
pus clouds are located between the Upper-Centaurus- 
Lupus and Upper-Sco subgroups of the Sco-Cen OB 
association, apparently lying on the edge of an Hi shell 
blown by the latter (and younger) of these (de Geus 
1992). Parts of the cloud complex show evidence of 
interaction with the expanding shell, and it has been 
suggested that star formation may have been triggered 
by both the current interaction and a previous shock 
wave from the older Upper-Centaurus-Lupus shell sev- 
eral Myr ago (Tachihara et al. 2001; Merin et al. 2008; 
Tothill et al. 2009). An implication of this is that the 
molecular matter - or at least some of it - pre-dated 
these shells. Along similar lines, the detailed work on 
the Ophiuchus region by de Geus (1992) stresses the 
interaction of the Upper-Sco subgroup with the much- 
studied p-Oph cloud. Here, molecular gas clearly pre- 
dates the current interaction, and is being stripped 
from the cloud to form streamers and tails. However, 
he also suggests the possibility of in-situ formation in 



the shell walls for some other parts of the complex. 
This is consistent with an alternative scenario sug- 
gested by Preibisch & Zinnecker (2007), in which feed- 
back flows formed both the present Lupus and Ophi- 
uchus complexes, and the Lupus clouds were formed 
sandwiched at the interface where the two shells meet. 
This scenario is attractive in that it does not require 
molecular gas to exist for an unrealistically long time 
prior to the onset of star formation. The authors also 
stress the role of multiple episodes of energy input, 
with clouds initially being formed during the wind 
bubble stage and then swept over to trigger star forma- 
tion during a later supernova phase. It is interesting to 
note that this is also a solution proposed by Inoue & 
Inutsuka (2009) to overcome the difficulty of achiev- 
ing high enough densities for molecular gas and star 
formation in a magnetically supported cloud. 

Another interesting Gould Belt complex is the Cepheus 
Flare. This region contains ~ 5 x 10 3 Mq of molecular 
gas that forms part of an expanding shell of R ~ 50 pc 
enclosing an old SNR (Grenier et al. 1989; Olano et al. 
2006). Kun (1998) note that the location of star form- 
ing sites at the edges of the CO clouds suggest external 
shock-triggering, which Olano et al. (2006) attribute to 
the same energetic event that created the shell - again, 
tentatively suggesting that the clouds pre-date the cur- 
rent interaction. Interestingly, this shell may be asso- 
ciated with the prominent Hi ring known as the North 
Celestial Pole Loop (Meyerdierks et al. 1991), which it- 
self contains the well-known molecular cirrus cloud, the 
Polaris Flare (Heithausen & Thaddeus 1990). The Po- 
laris Flare is a classic example of a translucent molec- 
ular cloud - diffuse (Meyerdierks & Heithausen 1996), 
gravitationally unbound (Heithausen & Thaddeus 1990) 
and non-star-forming (Andre et al. 2010). It has been 
described as 'the archetype of the initial phases of 
molecular cloud formation' (Miville-Deschenes et al. 
2010). In this context it is especially interesting to 
note that the Polaris Flare may be material swept up 
by the Cepheus Flare shell. 

Somewhat further afield, Heiles (1998) suggest a 
scenario in which star formation in the Orion and Gum 
regions (both at D ~ 500 pc) was triggered by a more 
distant supershell GSH 238+00+09; although it is un- 
clear whether molecular clouds were formed by the ac- 
tion of the shell or merely compressed by it. This sec- 
ondary star formation later blew the Orion-Eridanus 
(e.g. Heiles et al. 1999) and Gum Nebula (e.g. Yam- 
aguchi et al. 1999a) bubbles, both of which are them- 
selves associated with CO-bright gas, including many 
cometary globules - presumably remnants the parent 
clouds of the new young clusters, (see Bally et al. 1998; 
Yamaguchi et al. 1999a, and references within). 

5.3 The Milky Way Galaxy 

In the wider Milky Way, a large number of studies have 
documented molecular gas associated with supershells 
(Kundt & Mueller 1987; Koo & Heiles 1988; Heiles 
et al. 1996; Maciejewski et al. 1996; Normandeau et al. 
1996; Jung et al. 1996; Bally et al. 1998; Patel et al. 
1998; Rizzo & Arnal 1998; Yamaguchi et al. 1999a; 
Fukui et al. 1999; Kim & Koo 2000; Carpenter et al. 
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Figure 4: Subsection of the wall of the Galactic supershell GSH 287+04-17 integrated over the velocity range 
indicated in the top-left corner of the image (reproduced from Dawson et al. 2012a). The greyscale is Hi and the 
filled contours are 12 CO(J=1-0). Features labelled a indicate small (Mh 2 ~ 1O 2_3 M0) clumps of molecular gas 
offset towards the tips of atomic fingers that point in the direction of the shell centre. These are CO clouds that 
are likely being destroyed by their interaction with the shell. The feature labelled b is an example of a larger 
(Mh 2 ~ 1O 4 M0) molecular cloud that is well embedded in atomic material, and forms a coherent part of the main 
shell wall. This is a strong candidate for in-situ formation of molecular gas in the shell wall. This entire region is 
located at z ~ 200-300 pc above the Galactic midplane. 



2000; McClure-Griffiths et al. 2000; Matsunaga et al. 
2001; Moriguchi et al. 2002; Megeath et al. 2003; Moor 
& Kiss 2003; Dawson et al. 2008b, 2011b), including 
notable examples of multiple generations of sequential 
star formation triggered by large-scale stellar feedback 
(Patel et al. 1998; Oey et al. 2005). However, as in 
the case of the local ISM, there is often ambiguity sur- 
rounding the origins of the molecular material itself. 
Below we will discuss those objects for which the issue 
has been explicitly considered. 

The Perseus Chimney is a vertically elongated (~ 
110 x 230 pc) supershell in the Perseus Arm of the 
Outer Galaxy, which is likely in the process of blow- 
ing out of the Galactic Disk (Normandeau et al. 1996; 
Dennison et al. 1997; Basu et al. 1999). The shell is 
bordered on its lower edge by the bright W3/W4 Hn 
regions, and is associated with substantial quantities 
of molecular gas, particularly around its base in the 
Galactic Plane (Digel et al. 1996; Heyer & Terebey 
1998). Heyer et al. (1996) and Taylor et al. (1999) re- 
port on two cometary molecular clouds embedded in 
the chimney cavity, at altitudes of 20 pc and 50 pc and 
with molecular masses of ~ 5 x 10 3 and ~ 3 x 10 3 Af . 
These clouds show strong evidence of interaction with 
the stellar winds and UV radiation of the IC 1805 clus- 
ter at the chimney base, including streamers and tails 
of both Hi and CO pointing directing away from the 
ionising sources. They are both well interpreted as 
remnants of an original molecular complex that was 
destroyed in the formation of the chimney. For the 



large (~ 10 s Mq; Digel et al. 1996) reservoir of molec- 
ular gas around the base of the supershell, interpre- 
tations are less clear. Star formation along the edge 
of the shell in material bordering the W4 region im- 
plies compression and triggering, with the most natu- 
ral explanation being that this occurs in pre-existing 
molecular gas (Carpenter et al. 2000; Oey et al. 2005). 
However, whether all of this material represents rem- 
nants of an older GMC complex, or whether additional 
episodes of compression have supplemented the molec- 
ular gas reservoir in the past is not clear. Some hint 
that this might have been the case comes from the sug- 
gestion of Oey et al. (2005) that the IC 1805 cluster 
itself may represent a later generation of star forma- 
tion. In this picture its formation was triggered by 
an older progenitor cluster responsible for blowing an 
even older and larger supershell of which the current 
Perseus Chimney forms only the lower part (Reynolds 
et al. 2001). 

Several other studies have considered molecular cloud 
formation in a more quantitative way. Jung et al. 
(1996) report the association of a large mass of molec- 
ular gas (~ 1.1 x 10 6 Mq) with the Outer Galaxy 
supershell GS 234-2, which they interpret in the con- 
text of the gravitational fragmentation of the swept-up 
atomic medium based on the analytical expression of 
McCray & Kafatos (1987). Kim & Koo (2000) make a 
more careful analysis of two small (~ 10 3 Mq) molecu- 
lar clouds well-embedded within the Hi of the 'Galac- 
tic Worm' GW 46.4+5.5, which makes up the verti- 
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cal wall of a large (340 x 540 pc) and relatively local 
(D ~ 1.4 kpc) Hi shell. They also use the analytic ex- 
pression for gravitational fragmentation, and find that 
the timescale for its onset is approximately equal to 
the estimated kinematic age of the shell (t ~ 5 Myr), 
leading them to conclude that the molecular gas has 
likely formed from the swept-up ambient medium. In 
this context it is interesting that the molecular clouds 
(as measured from CO) do not appear to be globally 
gravitationally bound. While the molecular portion of 
the ISM presumably does not account for the entire 
mass of the 'fragment' in which it is found, it is worth 
noting that unbound clouds are not inconsistent with 
the predictions of the colliding flows picture of molec- 
ular cloud formation. 

5.3.1 GSH 287+04 17 and GSH 277+00+36: 
Detailed Case Studies and Quantita- 
tive Analysis 

The most detailed analysis of the origin and evolution 
of molecular gas in Milky Way supershells is given in 
a series of papers by Dawson et al. (2008a,b, 2011a,b). 
These authors present matched resolution, parsec-scale 
observations of Hi and CO in two Galactic supershells, 
GSH 287+04-17 and GSH 277+00+36, enabling de- 
tailed investigation of the relationship between the atomic 
and molecular ISM in shell walls. They find substan- 
tial quantities of co-moving molecular gas in the Hi 
shells, with rich substructure in both tracers, includ- 
ing molecular gas seen elongated along the inner edges 
of the atomic walls, embedded within atomic filaments 
and clouds, or taking the form of small CO clouds at 
the tips of tapering atomic 'fingers.' Figure 4 shows a 
section of the wall of GSH 287+04-17 that illustrates 
these features. Small CO clouds at the tips of Hi fin- 
gers have no substantial atomic envelopes and show no 
evidence for hidden 'dark' material (either optically 
thick Hi or CO-dark H2), implying that insufficient 
material exists to shield CO against photodissociation. 
Some also show shock-disrupted morphology, leading 
to the interpretation that these small clouds are pre- 
existing molecular gas currently undergoing dynami- 
cal disruption, gas stripping, and eventual dissociation 
due to their interaction with the shell. Survival life- 
times are roughly estimated as a few Myr, both from 
reference to the numerical results described in §4.1 and 
through estimates of the photodissociation timescale. 

Conversely, CO clouds well embedded within the 
main atomic shell walls are excellent candidates for 
newly formed clouds. The feature labelled 'b' in the 
figure is the strongest example of such an object. This 
cloud shows evidence for a substantial dark component 
identified from lOOpim IR excess, which comprises over 
50% of the total mass of the Hl-CO complex, and pro- 
vides sufficient material to shield CO molecules against 
UV dissociation, demonstrating that the cloud can sur- 
vive and continue to grow in its present form. The 
mean initial density required to sweep up the complex 
from the ambient medium is ~ 1-10 cm~ 3 (depend- 
ing on the assumed geometry), consistent with a re- 
alistic mixture of WNM and CNM. Comparing with 
the numerical simulations of §3.3, and assuming a flow 



speed equal to the present-day expansion velocity of 
v exp ~ 10 km s _1 , this implies that timescales of < 10 7 
yr are needed to form significant quantities of CO - 
consistent with the estimated age of the shell. In ad- 
dition, Wiinsch et al. (2012) argue that the mass spec- 
trum of molecular clumps in this region is consistent 
with the predictions of pressure-assisted gravitational 
instability in an expanding shell (PAGI; Wiinsch et al. 
2010), providing further constraints on the physics of 
the swept-up ISM. 

The net effect of the two shells on the volumes 
they occupy is estimated by comparing the molecular 
fraction, /h 2 = Mh 2 /(Mhi + Mh 2 ), in shell volumes 
- including the evacuated voids - to that in their lo- 
cal vicinities (essentially a proxy for the undisturbed 
medium) . Since /h 2 varies with location in the Galac- 
tic Disk, these 'background' regions are restricted care- 
fully to include only emission that is genuinely local to 
the shells. The results of this analysis reveal that /h 2 
in the supershell volumes is enhanced by a factor of 
~ 2-3 with respect to their local surroundings, imply- 
ing that as much as 50-70% of the molecular matter in 
their walls may have been formed as a direct result of 
stellar feedback. At present this analysis is restricted 
to two objects - both of which were selected in part be- 
cause of their association with molecular clouds - and 
no strong conclusions can be drawn about the Milky 
Way as a whole. Nevertheless, this is a compelling 
proof of concept, and some of the first quantitative 
evidence of molecular cloud formation in shell walls. 

5.3.2 Cold Hi in Shell Walls 

The presence of cold atomic medium in supershell walls 
is important, since it is a necessary precursor to molec- 
ular gas. Heiles (1982) report measurements of Hi ab- 
sorption spectra in the walls of several nearby super- 
shells, taken towards background continuum sources. 
They find that Hi in the shells is in the cold, neutral 
regime, with excitation temperatures ranging from 35- 
200 K - much colder than temperatures outside the 
shell walls. 

Dawson et al. (2011b) fit multiple Gaussian profiles 
to Hi emission in selected subsections the walls of the 
two shells discussed in §5.3.1, and use the measured 
linewidths to put a strict upper limit on the kinetic 
temperature of the gas of ~ 350 K. Assuming turbu- 
lence contributes approximately half of the observed 
linewidth (e.g. Heiles & Troland 2003), an estimate of 
Tk ~ 100 K is found for the fitted sections of shell wall. 
Together with an estimated density of ~ 10 cm~ 3 this 
implies that the atomic shell walls are dominated by 
cold gas with parameters close to the canonical values 
for the CNM (see also McClure-Griffiths et al. 2003). 

Hi self-absorption (HISA) provides a better probe 
of the morphology and physical properties of cold atomic 
clouds, although the requirement of bright background 
emission means that HISA is typically only observable 
at low Galactic latitudes. Moss et al. (2012) report the 
discovery of the very large supershell GSH 006-15+7, 
whose lower regions are seen as a striking cone-like ab- 
sorption feature against the Galactic Plane. They esti- 
mate optical depths and kinetic temperatures of r ~ 3 
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and Tk ~ 40 K, respetively, again indicating the pres- 
ence of copious quantities of cold, opaque atomic mate- 
rial. A similar and more extreme result is reported by 
Knee & Brunt (2001), who present evidence for an ex- 
tremely massive (~ 10 7 M Q ), extremely cold (~ 10-20 
K) atomic arc believed to belong to the distant outer 
galaxy shell GSH 139-03-69. 

Although little work on the thermal state of the 
atomic medium in supershells has been carried out be- 
yond these studies, they provide strong evidence that 
much of the ISM in supershells is in the form of cold 
atomic gas, consistent with them playing an important 
role in ISM cooling and cloud formation. 



0.5 

•I 0.4 

o 
to 

t 0.3 



% 0-2 



0.1 




♦ 

♦ 


♦ Dawson etal. (201 1a) 

■ Jung etal. (1996) 
m Patpl pt a f1 QQRl 

r cuo i ci ci I . y I jju| 

♦ Rizzo etal. (1998) 

■ Moor etal. (2003) 

♦ Moriguchi etal. (2003) 
x Yamaguchi etal. (1999) 


• 

• 

■ 


■ 

X 

• 



7 8 9 10 

Galactocentric Radius (kpc) 



5.3.3 High-altitude Molecular Material 

Another role of supershells may be in supporting a 
molecular 'thick disk' (Dame & Thaddeus 1994; Mal- 
hotra 1994a) by elevating molecular material well above 
its normal scale height (<r z ~ 60 pc at the solar circle 
Malhotra 1994b). Several studies have explicitly noted 
high-altitude molecular gas associated with feedback 
superstructures. Megeath et al. (2003) report the as- 
sociation of the large, high-altitude (z ~ 300 pc) star- 
forming molecular cloud complex NGC281 with an ex- 
panding supershell. Both of the shells described in 
detail in §5.3.1 also contain moderately sized (~ 10 3 - 
10 4 Mq) molecular clouds at heights of up to z ~ 450 
pc above the midplane (Dawson et al. 2011b). This in- 
cludes objects nominally identified both as pre-existing 
and newly formed clouds. 

Matsunaga et al. (2001) also report the detection 
of eight supershell candidates initially identified from 
CO alone, mostly seen either as holes in the CO distri- 
bution or as strings of discrete molecular clouds that 
trace arc-like shapes above the plane. Six of these 
objects contain molecular gas at z > 150 pc. Mod- 
erately elevated clouds are also seen by Kim & Koo 
(2000) and Moriguchi et al. (2002) in GW 46.4+5.5 
and the M16/M17 supershell, respectively. Finally, 
Heiles et al. (1996) note the presence of 'high-z CO' as- 
sociated with two Galactic Worms, GW 30.5-2.5 and 
GW 49.1-1.4, though they do not quote z itself. These 
results suggest that the molecular thick disk is sup- 
ported at least in part by both the formation and the 
displacement of molecular gas due to discrete large- 
scale stellar feedback events. 

5.3.4 Remarks on Molecular Fraction and 
Galactocentric Radius 

The picture of supershell-driven molecular cloud for- 
mation outlined in this review applies primarily to a 
regime in which the ambient medium contains a large 
atomic component. Mean densities must also be high 
enough that feedback flows are able to accumulate 
sufficient quantities of material for gas cooling and 
molecule formation. These conditions are applicable 
within a few kpc of the solar circle, but are unlikely 
to apply to the extreme inner or outer regions of the 
disk. The relevance of feedback-triggered molecular 
cloud formation is therefore likely to be a function of 
Galactocentric radius. 



Figure 5: The relationship between Galactocentric ra- 
dius and molecular fraction in Milky Way supershells. 
The molecular fraction, /h 2 = Mu 2 / (Mm + Mu 2 ) is 
calculated for all supershells in the literature for which 
estimates of both atomic and molecular mass are avail- 
able for the full shell. 



While quantitative measurements of molecular cloud 
formation in supershells are still too few to provide any 
observational insight into its variation with location, it 
is worth noting that available estimates of the Hi and 
H2 masses of Galactic supershells do demonstrate the 
expected propensity for increasing molecular fraction 
with decreasing Galactocentric radius (see figure 5). 
This also has implications for the initial detection and 
statistical study of populations of supershells, particu- 
larly when detections are based on a single tracer. The 
eight supershell candidates discovered in CO by Mat- 
sunaga et al. (2001) are all inner Galaxy objects, lo- 
cated between R ~ 5.5-7.5 kpc (rescaled for Ro = 8.0 
kpc). Conversely, Hi-based detections may well be bi- 
assed against inner Galaxy shells, as decreasing ambi- 
ent atomic gas fraction and increasing confusion ren- 
ders coherent Hi structures more difficult to detect. 

5.4 The Large Magellanic Cloud 

The Large Magellanic Cloud (LMC) is the nearest star- 
forming galaxy to the Milky Way (D ~ 50 kpc), and 
an excellent laboratory for studying the interaction be- 
tween stars and the ISM. Unlike in the Milky Way, 
line of sight confusion is not a significant problem, and 
the reliable identification of supershells is considerably 
simplified, enabling a more statistical approach to the 
question of molecular cloud formation. A large num- 
ber of shells and supershells have been identified (e.g. 
Davies et al. 1976; Meaburn 1980; Chu & Mac Low 
1990; Kim et al. 1999), in a variety of observational 
tracers, and several studies have explicitly considered 
their relationship to the molecular ISM. 

Yamaguchi et al. (2001b) make a statistical study 
of the effects of optically identified supergiant shells 
(SGSs; defined as those whose radii are larger than the 
disk scale height) on the formation of stars and molec- 
ular clouds. They find that both the number density 
and mass density of 12 CO(J=1-0) molecular clouds is 
enhanced by a factor of 1.5-2 at the edges of SGSs, 
which they interpret as evidence of cloud formation. 
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Figure 6: Supergiant shells and shell complexes overlaid on integrated intensity map of the LMC (reproduced 
from Dawson et al. 2012b). The greyscale image is Hi (Kim et al. 2003), and blue contours are 12 CO(J=1-0) 
(Fukui et al. 2008), processed as described in Dawson et al. (2012b), and integrated over the full velocity range of 
the LMC disk. The solid blue line marks the boundary of the region observed in CO. Dark pink lines trace the 
inner rims of the shell complexes and purple lines mark their outer boundaries (delineating the outer edges of the 
dense shells). Dotted blue lines enclose the region known as the southeastern Hi overdensity. 



Book et al. (2009) examine a subset of this SGS pop- 
ulation in detail, and argue that both the compression 
of pre-existing dense clouds and the formation of new 
molecular gas are likely occurring in the supergiant 
shells LMC 1, 4, 5, and 6. The interacting shells LMC 
4 and LMC 5 contain particularly convincing examples 
of both processes: A large, dense ridge of molecular 
material compressed between them forms a striking ex- 
ample of the classic picture of efficient cloud formation 
at the interface of two shells (e.g. Hartmann et al. 2001; 
Ntormousi et al. 2011), while a smaller cometary CO 
cloud with an Hi tail exhibits the classic morphological 
signatures of the interaction of a shock front with pre- 
existing dense gas (see also Yamaguchi et al. 2001a). 
Hi absorption spectra towards background continuum 
sources also suggest that efficient cooling has driven 
the production of large quantities of cold atomic gas 
in the walls of LMC 4 (Marx-Zimmer et al. 2000). This 



is consistent with the behaviour seen in Milky Way su- 
pershells. 

Dawson et al. (2012b) combine CO data with Hi 
synthesis images (Kim et al. 2003) to make the first 
quantitative measurement of feedback-driven cloud for- 
mation in an entire galactic system. As described in 
§5.3.1 for Milky Way shells, their method compares 
the molecular fraction, /h 2 = Mn 2 /(Mm + Mn 2 ), in 
SGS volumes with that elsewhere in the disk, in order 
to assess whether the passage of an SGS through the 
ISM is associated with a net increase in the production 
of molecular gas. Figure 6 shows an Hi and CO map 
of the LMC disk, overlaid with the outlines of the SGS 
complexes used in this analysis. They find that /h 2 in 
the aggregate volume occupied by all SGSs is identical 
to the rest of the LMC disk, suggesting that supergiant 
shells are not a dominant driver of molecular cloud for- 
mation on galactic scales. Indeed, the global structure 
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of the LMC disk is found to be a better determinant of 
where the highest molecular fractions are found. How- 
ever, the majority of objects (~ 70% by mass) are 
more molecular than their local surroundings, imply- 
ing that the presence of a supergiant shell does on av- 
erage have a small positive effect on the molecular gas 
fraction. This analysis is used to place a lower limit on 
the total fraction of molecular cloud formation in the 
LMC driven by large-scale stellar feedback, which is 
estimated to be > 4-11% of the total molecular mass 
of the galaxy. 

The importance of stellar feedback as a driver of 
molecular cloud formation is expected to be dependent 
on galaxy type. As a dwarf irregular, the Large Magel- 
lanic Cloud (LMC) is expected to be more susceptible 
to feedback-triggered molecular cloud formation than 
early-type spirals such as the Milky Way (Elmegreen 
et al. 2002, see also §3.2). With low shear (Weidner 
et al. 2010), a large Hi scale height (Brinks et al. 2002) 
and weak spiral structure, shells in the LMC are able to 
grow larger before they are deformed, expand further 
before vertical blowout and depressurisation, and not 
suffer the spiral arm disruption that will affect their 
Milky Way counterparts. Similarly, the lack of strong 
spiral structure, combined with a relatively weak disk 
gravitational potential, suggests that galaxy-scale self 
gravity and the accumulation of the ISM in spiral arms 
may play a less important role in dense gas formation 
than in grand-design spirals. 

6 Summary, Discussion and Fu- 
ture Directions 

The accumulation, compression and cooling of the ISM 
in supershells formed by OB cluster feedback can drive 
the production of star-forming molecular clouds. Ana- 
lytical and numerical models provide a strong theoret- 
ical basis for this process, and observations leave little 
doubt that molecular cloud formation via large-scale 
feedback is occurring in both the Milky Way and the 
LMC. 

Theoretically, many details remain to be hammered 
out, including the ability of magnetic fields to prevent 
compressed material from reaching sufficient densities 
for molecule formation, whether multiple episodes of 
compression are necessary, and whether the initial am- 
bient density must be higher than canonical WNM val- 
ues in order for flow-driven cloud formation to occur, 
(e.g. Dobbs et al. 2012; Inoue & Inutsuka 2012). Ob- 
servationally, there are many convincing candidates for 
molecular clouds formed in supershell walls, and am- 
ple evidence that the atomic ISM in these shell walls 
tends to be cold, demonstrating that the cooling of 
swept-up material occurs readily. Enhanced molecular 
gas fractions in the volumes occupied by two Galactic 
supershells have also provided one of the first quantita- 
tive measurements of feedback-driven cloud formation 
in the Milky Way (Dawson et al. 2011a). However, 
the question of how readily the process occurs remains 
open, as does its relative importance in the galactic 
context. The first direct measurement of the contribu- 
tion of stellar feedback to molecular cloud formation 



on galactic scales has now been made in the LMC, 
where > 4-11% of the total molecular gas content of 
the galaxy is estimated to have been formed as a direct 
result of large-scale feedback (Dawson et al. 2012b). 
However, no such estimates exist for the Milky Way or 
other nearby galaxies. 

Thinking globally, large-scale gravitational insta- 
bilities play a major role in the initial formation of 
dense gas in galaxy disks (see e.g. Elmegreen 2002; 
Mac Low & Glover 2012, and references within), and 
may be the primary driver of molecular gas forma- 
tion in most systems. Spiral arms are also clearly an 
important means of accumulating material, although 
whether they drive a significant amount of molecular 
gas formation, or whether their primary influence is 
to agglomerate existing GMCs into larger complexes is 
not yet clear (Koda et al. 2009). A picture is suggested 
in which large-scale feedback plays a secondary role - 
initiating further episodes of compression and addi- 
tional cloud formation in an inhomogeneous medium 
that already contains significant amounts of dense struc- 
ture. This picture is tentatively consistent with the 
results presented in this review. However, direct mea- 
surements of the contribution of stellar feedback to 
molecular cloud formation rates are strongly desirable. 

Observational work that explicitly deals with feedback- 
triggered cloud formation has generally taken a morphology- 
based approach, which relies on first identifying su- 
pershell structures and then examining the associated 
molecular ISM. This strategy requires that the neu- 
tral ISM be resolved to scales of ~ 10-100 pc, which 
has historically been challenging for external galax- 
ies, particularly in molecular line observations. With 
the advent of millimetre and sub-millimetre interfer- 
ometers such as the Combined Array for Research in 
Millimeter Astronomy (CARMA) and the Atacama 
Large Millimeter/submillimeter Array (ALMA), high- 
resolution, high-sensitivity observations of molecular 
tracers are now becoming possible in local group galax- 
ies. This opens up the exciting possibility of extend- 
ing the study of feedback-triggered molecular cloud 
formation to large samples of galaxies with a range 
of different properties. Such work will be essential in 
constructing a consistent model of the impact of stel- 
lar feedback on ISM evolution across a broad range of 
galaxy types. 

Gaseous disks contain large amounts of stochastic 
structure, and the reliable identification of signatures 
of stellar feedback - particularly in the neutral ISM 
alone - is a source of uncertainty in observational stud- 
ies. A preferable strategy would be a means of explor- 
ing feedback-triggered cloud formation that does not 
first rely on the morphological identification of shell- 
like structures. A promising approach is suggested by 
the work of Dobbs et al. (2012), who identify character- 
istic distributions in the eigenvalues of the local rate- 
of-strain tensor in model disks dominated by different 
astrophysical drivers of converging flows. This sug- 
gests the possibility of identifying the primary drivers 
of the flows that assemble dense clouds without re- 
course to morphological arguments, and hence quan- 
tifying the importance of stellar feedback in a robust 
and statistical sense. While it remains to be seen if this 



16 



Publications of the Astronomical Society of Australia 



method can be usefully applied to real astronomical 
datasets (which cannot directly provide information 
on the 3D velocity field of the ISM), such direct mar- 
riages between simulations and observations are likely 
to provide powerful tools for the interpretation of fu- 
ture high-resolution data. This underscores the im- 
portance of using model galaxies - in which the input 
parameters are well constrained - to develop observa- 
tional diagnostics for the ISM in real systems. 
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